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Abstract 


An  investigation  was  conducted  to  determine  the 
effects  of  spring  and  inertia  devices  on  the  longi- 
tudinal static  and  dynamic  stability  of  aircraft. 
Theoretical  studies  and  flight  tests  of  a modified 
Navion  airplane  were  completed.  It  was  concluded 
that  introduction  of  bobweights  in  the  elevator  con- 
trol system  consistently  improved  the  phugoid  damp- 
ing and  downsprings  had  a favorable  effect  when  lim- 
ited to  a light  or  moderate  application. 
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I.  SUMMARY 


Theoretical  studies  and  flight  test  measurements  of  the  effects  of  bob- 
weights  and  downsp rings  on  static  and  dynamic  stability  have  been  completed. 

An  analytical  analysis  shows  the  springy  tab  is  equivalent  to  the  downspring. 
The  test  article  was  a modified  Navion  type  airplane  considered  in  cruising 
flight.  Theoretical  results  based  on  estimated  characteristics  were  initially 
obtained  for  complete  ranges  of  control  device  forces  and  center  of  gravity 
displacements  in  order  to  establish  qualitatively  the  significant  trends  to 
be  verified  by  arctual  measurement.  A flight  region  selected  on  the  basis  of 
the  analytical  computations  was  then  explored  in  a series  of  dynamic  tests. 
Static  tests  were  conducted  in  conjunction  with  this  phase  to  determine  the 
reference  static  stability  for  each  dynamic  test  condition.  Random  checks 
of  the  dynamic  runs  were  also  made  with  an  analog  computer  utilizing  aerody- 
namic data  extracted  from  the  static  flight  results. 

The  analytical  phase  indicated  that  the  effects  of  force  devices  on  the 
longitudinal  dynamics  could  be  expressed  as  a function  of  their  control  of 
just  two  principal  static  parameters  — the  stick-free  maneuvering  margin  and 
the  difference  between  the  stick-free  maneuver  and  neutral  points.  The  first 
is  regulated  by  bobweight  and  center  of  gravity  location  and  the  second  is 
regulated  by  downspring  or  springy  tab  solely.  The  increment  /V^~  A/J  is 
referred  to  as  the  static  damping  parameter  since  this  term  is  dependent  upon 
the  basic  damping,  in  pitch  for  airplanes  without  auxiliary  control  devices 
and  in  the  subject  investigation  it  was  found  to  have  a strong  influence  on 
the  danping  characteristics.  It  was  established  that  the  phugoid  mode  was  of 
predominant  importance  in  the  static  stability  region  under  study.  All  per- 
tinent phugoid  trends  were  substantiated  by  the  flight  data  and  a reasonable 
agreement  was  shown  with  regard  to  the  magnitude  of  the  danping  and  period. 

In  addition  qualitative  comparisons  of  the  analog  solutions  and  the  flight 
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results  were  satisfactory. 

In  general,  it  was  found  that  increasing  the  stick-free  static  stability 
at.  rearward  displaced  centers  of  gravity  by  introducing  bobweight  forces  had 
consistently  favorable  effects  on  the  phugoid.  The  mode  changed  from  a pure 
divergence  to  a damped  oscillation  at  low  values  of  downspring  force  (positive 
damping  parameter)  and  from  a pure  divergence  to  an  undamped  oscillation  at 
high  values  of  downspring  force.  A moderate  application  of  downspring  to  re- 
store level-flight  static  stability  improved  the  phugoid  damping  from  a pure 
divergence  to  a damped  oscillation  when  the  spring  forces  were  limited  to  main- 
tain a positive  damping  parameter.  Further  increases  in  downspring  force  suf- 
ficient to  cause  the  damping  parameter  to  become  negative  resulted  in  rapidly 
divergent  oscillations.  This  effect  was  particularly  pronounced  at  low  posi- 
tive maneuver  margins.  In  all  cases  a damped  phugoid  was  obtained  only  for 
a positive  damping  parameter  in  conjunction  with  a moderate  maneuver  margin. 

Practical  limits  for  application  of  bobweights  and  downsprings  were 
found  to  depend  upon  the  basic  elevator  and  tab  power  and  thereby  imposed  re- 
strictions on  the  center  of  gravity  displacement  which  could  be  corrected  by 
these  control  modifications.  Other  considerations  made  evident  in  the  test 
program  were  the  possibility  of  adverse  effects  due  to  the  increases  in  ele- 
vator inertia  and  ground-handling  stick  forces  that  occur  when  auxiliary  con- 
trol devices  of  considerable  strength  are  used. 

Most  of  the  flight  testing  was  conducted  with  the  center  of  gravity 
located  aft  of  the  stick-fixed  neutral  point  so  that  the  airplane  had  an  un- 
stable elevator  deflection  gradient.  It  is  of  interest  to  note  that  through- 
out the  program  the  pilot  had  no  impression  in  the  handling  qualities  of  this 
static  instability.  The  concept  that  only  the  stick-free  stability  has  im- 
portance in  level-flight  operation  would  require  further  comprehensive  test- 
ing for  confirmation. 
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II.  INTRODUCTION 


Modern  trends  in  the  design  of  high  performance  aircraft  have  made  it 
increasingly  difficult  to  achieve  center  of  gravity  positions  located  suffi- 
ciently forward  to  meet  requirements  for  satisfactory  static  stability  force 
gradients.  Therefore,  artificial  means  of  improving  static  stability,  although 
used  fairly  extensively  in  the  past,  have  become  more  prevalent.  The  stability 
change  is  accomplished  by  introducing  constant  forces  in  the  longitudinal  con- 
trol system  with  devices  such  as  bobwieghts,  downsprings  and  springy  tabs.  The 
usable  center  of  gravity  range  is  generally  limited  by  the  maximum  allowable 
maneuvering  force  gradient  and  the  minimum  allowable  level-flight  force  grad- 
ient. Through  the  application  of  downsprings  or  springy  tabs,  the  usable 
center  of  gravity  range  can  be  increased  in  magnitude  and  with  bobweights  the 
absolute  values  of  each  gradient  can  be  increased  or  restored.  In  the  past, 
effects  of  these  devices  on  the  longitudinal  dynamic  stability  were  not  con- 
sidered particularly  important  since  generally  the  basic  airplane  character- 
istics were  such  that  the  short  period  mode  was  heavily  damped  and  the  phugoid 
period  was  long  enough  to  make  its  damping  relatively  inconsequential.  Current 
design  trends  however,  have  also  led  to  numerous  instances  of  unsatisfactory 
phugoid  damping  so  that  now  a new  need  exists  for.  information  on  how  the  lon- 
gitudinal dynamics  are  modified  by  auxiliary  control  devices.  Since  relatively 
little  was  known  about  these  effects  this  program  was  undertaken  in  an  effort 
to  provide  actual  test  data. 

It  was  thought  feasible  to  use  a small  light  airplane  to  establish  gen- 
eralized results  in  a qualitative  form.  Although  quantitative  effects  of 
bobweights  and  downsprings  depend  upon  the  particular  characteristics  of  the 
airplane  under  consideration,  trends  could  be  determined  which  would  be  of 
value  in  application  to  subsonic  non-boost  control  aircraft  in  general.  Re- 
sults of  such  an  investigation  would  be  useful  in  the  design  of  force  additive 
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devices  for  the  satisfaction  of  static  requirements  as  well  as  to  afford  some 
means  of  correcting  inherent  dynamic  difficulties  encountered  without  force 
installations.  In  addition,  as  phugoid  requirements  for  satisfactory  flying 
qualities  are  not  as  yet  established,  information  on  the  practical  limits  for 
marginal  dynamic  characteristics  is  of  value. 
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III.  DESCRIPTION  OF  AIRPLANE 


The  flight  research  article,  a Navion  airplane,  is  a single-engine,  low- 
wing  monoplane  provided  with  a retractable  tricycle  landing  gear  and  an  ad- 
justable pitch  propeller.  A general  view  of  the  test  airplane  is  showi  in 
Figure  1 and  the  physical  characteristics  are  given  in  Appendix  I.  The  basic 
elevator  configuration  consists  of  a full-span  surface  without  aerodynamic 
balance.  Wheel-type  dual  controls  are  included  and  the  elevators  have  adjust- 
able trim  tabs  on  each  panel. 

In  order  to  minimize  the  destabilizing  center  of  gravity  displacements 
required  for  test  purposes  the  horizontal  tail  area  was  reduced.  This  was 
accomplished  by  decreasing  the  tail  span  inches  to  obtain  an  18$  reduction 
in  total  area.  The  modified  horizontal  tail  is  shown  in  comparison  with  the 
original  configuration  in  Figure  2.  After  modification,  the  elevator  was 
statically  unbalanced  with  a trailing  edge  down  hinge  moment  equivalent  to  a 
push  stick  force  of  approximately  8 lbs.  Adjustments  were  made  to  the  sta- 
bilizer incidence  setting  to  compensate  for  trim  changes  due  to  the  horizontal 
tail  redesign.  The  elevator  trim  tab  area  vra.s  increased  by  a trailing  edge 
extension  to  enable  trim  of  heavy'  force  devices  and  the  full  down  elevator 
deflection  was  increased  by  a few  degrees  to  improve  the  control  effective- 
ness for  dynamic  testing.  The  friction  level  in  the  longitudinal  control 
system  was  reduced  to  approximately  1 lb.  and  the  rear  seats  were  removed 
to  allow  installation  of  the  instrumentation. 
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IV.  WEIGHT  AND  BALANCE 


After  installation  of  all  instrumentation  and  modification  of  the  tail 
configuration  the  Navion  gross  weight  was  2680  lbs.  with  full  fuel,  oil  and 
a crew  of  two.  The  center  of  gravity  location  for  this  loading  was  at  27% 

MAC. 

It  was  estimated  that  center  of  gravity  locations  approaching  b0%  MAC 
would  be  desirable  for  the  test  program.  To  achieve  these  aft  positions  it 
was  necessary  to  install  fixed  ballast  in  the  fuselage  under  the  dorsal  fin 
as  well  as  transfer  movable  ballast  within  the  cabin  f ollowing  take-off. 

Fixed  ballast  was  increased  in  increments  of  20  lbs.  to  a maximum  of  80  lbs. 
and  movable  ballast  weighing  12$  lbs.  was  displaced  in  varying  increments 
through  a distance  of  8 ft.  during  flight  to  attain  center  of  gravity  loca- 
tions between  27%  and  39% • Thus,  the  maximum  take-off  weights  were  approxi- 
mately 2900  lbs.  in  the  course  of  the  program.  Center  of  gravity7  locations 
aft  of  39%  were  not  considered  reasonable  for  the  static  tests  based  on  con- 
siderations of  overload  gross  weights,  ground  handling,  structural  limitations 
and  ease  of  ballast  transfer.  Dynamic  tests  were  restricted  to  center  of 
gravity  locations  forward  of  3&%  MAC  due  to  a deficiency  of  down  elevator  de- 
flection reqiired  for  recovering  from  severe  nose-up  attitudes.  Records  were 

kept  on  balance  changes  in  order  to . accurately  compute  take-off  weights  and 

v 

center  of  gravity  locations  for  each  flight.  Effects  on  the  balance  of  fuel 
consumption  during  flight  were  found  to  be  negligible. 
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V,  BOBWEIGHT  AND  DOWN  SPRING  SYSTEMS 


A rigid  arm  was  fixed  normal  to  the  control  column  so  that  through  a 
reverse  leverage  system  a positive  (tending  to  depress  the  elevator)  bohweight 
moment  could  be  applied  and  by  attaching  a preloaded  spring  from  the  canopy 
to  this  member  a positive  downspring  moment  could  be  introduced  as  well.  The 
test  installation  is  shown  in  Figure  3*  Before  the  test  phase  was  begun  the 
elevator  was  statically  balanced  by  attaching  a weight  directly  to  the  con- 
trol column  arm  to  correct  the  trailing  edge  heaviness  condition.  Both  the 
bobweight  and  downspring  systems  were  designed  so  that  the  moment  exerted  on 
the  elevator  was  essentially  constant  in  the  test  range  of  deflections. 

Positive  bobweights  and  downspring  moments,  and  combinations  thereof, 
were  applied  progressively  in  equal  increments,  each  equivalent  to  a stick 
force  of  8 lbs.  This  was  done  using  a series  of  similar  weights  and  a series 
of  tension  springs  of  equal  preload.  The  magnitude  of  incremental  force  was 
selected  as  that  affording  a static  stability  variation  consistent  with  con- 
venient changes  in  airplane  balance. 
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VI.  INSTRUMENTATION 


Special  test  instrumentation  was  installed  in  the  Navion  to  accurately 
measure  airspeed,  control  force,  normal  acceleration  and  elevator  deflection. 
For  the  static  tests  readings  were  visually  recorded  and  for  dynamic  testing 
a photorecorder  was  used  including  a 16  mm  gun  camera  to  film  an  airspeed 
indicator  and  a stop  watch. 

Airspeed  data  were  obtained  with  a pitot-static  head  located  on  a wing 
boom.  Calibrations  indicated  the  corrections  for  instrument  and  position 
error  were  negligible.  The  control  force  applied  to  the  wheel  was  measured 
with  a strain  gage  system  to  give  current  readings  which  were  proportional  to 
wheel  force  and  accurate  to  0.5  lb.  The  accelerometer  calibrated  to  read 
directly  in  "g"  units  consisted  of  a freely  suspended  spring  and  mass  en- 
closed in  a glass  tube  and  was  accurate  to  .01  "g" . Elevator  .positions  were 
measured  to  within  0.1  degree  with  a dual  autosyn  system  made  up  of  a master 
unit  attached  through  a gear  train  to  the  elevator  torque  tube  which  drove  a 
repeater  indicator  in  the  cabin.  Stick  force  and  elevator  deflection  c ali- 
brations  were  kept  up  to  date  throughout  the  test  phase  to  minimize  error 
due  to  sensitivity  changes  and  drift. 

\ 


WADC  TR  53-350 


a 


VII.  THEORETICAL  ANALYSIS 

A.  Static  Stability 

The  longitudinal  stick-free  stability  of  aircraft  exhibiting  excessively 
low  control  force  gradients  may  be  improved  by  introducing  artificial  forces 
which  are  independent  of  airspeed.  This  is  accomplished  by  placing  devices 
such  as  bobweights  and/or  downsprings  in  the  longitudinal  control  system  to 
produce  constant  mechanical  hinge  moments  tending  to  increase  the  positive 
floating  angle  of  the  elevator.  Under  unaccelerated  flight  conditions  the 
stability  changes  resulting  from  each  device  are  equivalent.  When  the  air- 
plane develops  normal  accelerations  in  maneuvering  flight  the  bobweight  mo- 
ment, being  subject  to  inertial  effects,  increases  as  a linear  function  of 
the  acceleration  whereas  the  downspring  moment  remains  unchanged  from  the 
level-flight  condition.  The  static  stability  effects  of  bobweights  and  down- 
springs  are  readily  calculated  for  both  level- flight  and  maneuvering  flight 
operation. 

The  change  in  elevator  floating  angle  caused  by  the  control  devices  is 
given  by  / \ 

'***  Q (1) 

where  the  hinge  moment  coefficient  developed  mechanically  is 

— ,/>  s 

US  Ce 

The  airplane  pitching  moment  coefficient  corresponding  to  the  change  in 
floating  single  is 


(2) 


Cn(Ab). 


Combining  these  relations  yields 


far  i = fay)™ 

( ^ Cs  Q6 


(3) 


(K 
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Stability  contributions,  expressed  as  stick- free  neutral  point  displacements, 
are  now  evaluated  by  finding  the  rate  of  change  of  pitching  moment  with  lift 
coefficient.  The  lift  coefficient  is  defined  in  general  form  by 


(5) 


In  unaccelerated  flight  the  load  factor  is  maintained  constant  at 
and  the  airspeed  considered  as  a variable.  Consequently,  the  pitching  moment 
in  terms  of  lift  coefficient  is 


% Ce  Qs 


(6) 


and  taking  the  derivative  with  respect  to  gives 

= hL&S 

C^=/  W/S  J3e  Ce  Q6 


(7) 


This  expression  then  represents  the  shift  in  stick-free  neutral  point  pro- 
duced by  either  a bobweight  or  downspring  hinge  moment. 

In  maneuvering  flight  the  airspeed  is  held  fixed  at  some  trim  value  and 
the  angle  of  attack  or  normal  load  factor  is  varied.  Referring  to  the  general 
moment  equation  it  can  be  seen  that  for  the  downspring  there  is  no  change  in 
pitching  moment  for  this  condition  and  therefore  the  downspring  has  no  effect 
on  maneuvering  stability.  In  the  case  of  the  bobweight,  the  change  in  moment 
coefficient  is 


3 ^ Ce  C*6 


(8) 


and  the  derivative  is 


/dC»\  Qs±  „ 


The  rate  of  change  of  bobweight  hinge  moment  with  lift  coefficient  can  be 
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expressed  as 


^ f MS?  x <^22 

c2  Q ot-r?  c2C<_  ^ 


where  {///*? )^a)  _ /V,,^)  (ll) 

and  -5c-  = Tfe  <12> 

Substituting  to  obtain  the  final  form, 

/at-Gr,  V = 22. (13) 

^^ceQs 

This  term  then  represents  the  shift  in  stick-free  maneuver  point  due  to  the 
bobweight  and  is  seen  to  be  identical  to  the  shift  in  neutral  point  obtained 
with  either  a bobweight  or  downspring  of  equal  strength.  Therefore  it  may 
be  stated  for  future  reference,  the  bobweight  does  not  alter  the  basic  dif- 
ference between  the  neutral  and  maneuver  point  locations  whereas  the  down- 
spring does. 


Another  auxiliary  control  device  which  has  the  very  same  stability 
effect  as  the  downspring  is  the  springy  tab.  The  only  distinction  between 
the  two  is  the  fact  that  the  downspring  produces  a stick  force  noticeable  in 
ground  handling  and  the  springy  tab  has  no  effect  at  zero  airspeed.  This 
mechanism  was  not  evaluated  in  the  reported  study  since  it  was  demonstrated 
that  the  dynamics  depend  on  the  static  stability  changes  produced  by  the  de- 
vices. On  this  basis  if  the  springy  tab  is  properly  designed  it  should  have 
the  same  influence  on  the  dynamic  characteristics  as  the  downspring  when 
producing  the  same  change  in  static  stability.  The  static  equivalence  be- 
tween the  springy  tab  and  the  downspring  may  be  s hown  analytically. 

The  springy  tab  is  a free  elevator  tab  in  combination  with  a spring 
which  exerts  a constant  moment  about  the  tab  hinge  line.  At  zero  airspeed 


WADC  TR  55-350 


11 


the  tab  is  at  its  full  up  deflection  and  this  position  is  maintained  until  the 
the  airspeed  increases  to  a particular  value  where  the  spring  torque  is  bal- 
anced. For  speeds  in  excess  of  this  value  the  tab  will  deflect  downward  to 
some  equilibrium  floating  angle.  The  hinge  moment  about  the  tab  hinge  line 
is  given  by 

ML  - * + 4*  f-  ^ c*  a-) 

where  AZ  is  the  spring  torque  in  lb. -ft.  and  C,  is  the  rate  of  change 
of  tab  hinge  moment  with  tab  deflection.  Setting  this  expression  equal  to 
zero  and  solving  for  ^ yields  the  tab  floating  angle f 

The  elevator  hinge  moment  due  to  the  tab  deflection  is 

/AS?  ~ 

where  is  the  rate  of  change  of  elevator  hinge  moment  with  tab  deflec- 

tion.  Substituting  for  the  tab  deflection  yields 

/toy  = - — 7:  (17) 

This  expression  which  is  independent  of  airspeed  shows  that  the  elevator 
hinge  moment  produced  by  the  springy  tab  is  constant  just  as  is  that  of  the 
downspring.  And  when  the  springy  tab  installation  is  designed  so  as  not  to 
introduce  weight  moments  about  the  elevator  hinge  line  it  has  no  effect  on 
the  maneuvering  stability  as  in  the  case  of  the  downspring. 

B.  Dynamic  Stability 
1 . Exact  Solut ion 

The  general  longitudinal  equations  of  motion  of  an  airplane  with  the 
elevator  free  to  rotate  are  derived  in  reference  1.  In  non-dimensional  form 


(15) 

06) 


WADC  TR  53-350 


12 


the  four  linear  differential  equations  are 


Drag 

(Cj>  + c/)ic  - jt  (CL  ~ ClJ  »<.  & = o (18) 

Lift 

C,  ^ O ~ 0/0  = O 0-9) 

Pitching  Moment 

^ (C»  % c/)  6 = O 

Elevator  Hinge  Moment 

C\  *■  rfCi.-  U;  CuJo'A  +M'QJc/&i 

~ ( ^■('Qs  ^ c/^J  6 = O 

All  the  relations  are  referred  to  a time  base  of  <v/£-  and  the  operator  is 
defined  as 

The  qualifying  assumptions  involved  in  the  development  of  these  expressions 
are  that  the  motion  generated  xn  response  to  some  initial  disturbance  or  con- 
trol input  is  small  in  amplitude,  motions  depend  on  displacement  and  velocity 
dependent  derivatives  only  and  not  acceleration  terms  and  also  that  these 
partial  derivatives  or  aerodynamic  characteristics  are  constants  for  a given 
airplane  thus  making  the  differential  equations  linear.  These  are  major  con- 
siderations in  the  analytical  treatment  of  the  flight  dynamics  and  the  extent 
to  which  these  conditions  apply  determines  the  degree  of  agreement  with  the 
actual  airplane  motion. 

Appropriate  modifications  to  the  general  equations  are  first  made  for 


the  subject  study.  The  pitching  moment  due  to  the  rate  of  elevator  deflec- 
tion,  L_>»yg  > was  estimated  for  the  Navion  and  found  to  be  negligible.  The 
change  in  pitching  moment  with  flight  velocity,  , may  be  disregarded  at 
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this  point  since  compressibility  effects  are  not  involved  and  power  contri- 
butions to  the  pitching  moment  are  usually  small.  At  a later  stage  in  the 
analysis  it  is  ^shown  how  possible  power  effects  can  be  taken  into  account. 


Velocity  dependent  hinge  moment  terms  encompassed  by  (Zy  would  normally  be 
neglected  for  the  same  reasons.  However,  in  this  study  the  control  devices 
produce  an  effective  (2 a • The  equations  of  motion  as  given  are  for  no  in- 
puts;  that  is  they  represent  motion  in  return  to  trim  following  an  instanta- 
neous disturbance.  Hinge  moments  from  force  devices  constitute  a steady 
elevator  input  as  a right-hand  term  in  the  elevator  equation.  The  hinge  mo- 
ment derivative  is  then  a function  of  the  velocity  and  can  be  derived  con- 
sistent with  the  general  equations  as  follows: 

y v*.. % (22) 

Finding  the  partial  derivative  with  respect  to  velocity 


^ - q,  f ^ 


c 


(23) 


Dividing  by  (2^  an(^  multiplying  by  J/  yields  AL  <?t  as  the 

right-hand  term.  Transposing  in  the  elevator  equation,  it  is  apparent  that 
C = Ca  • After 


'k. 

become 

Drag 


making  the  preceding  changes  the  complete  equations 


*~c/ ) U £ v-  Z e 


--  O 


Lift 


(2k) 

(25) 


o u c.^c/y  ~ c/s  -c 

Pitching  Moment 

o' ) < + (C^  c/~  />(/*)&  <3(26) 

Elevator  Hinge  Moment 

A (6,  - Cj,h)c/j*  */(/i  * 

-(A*  A y-/:/ s y-  Ct^j-  4 c/J)d  = o (27) 


A simultaneous  solution  of  the  revised  equations  follows  the  standard 
procedure  given  in  reference  1 and  yields  a sixth  order  characteristic  equa- 
tion from  which  the  period  and  damping  of  each  of  the  three  modes  of  motion 
can  be  determined.  Usually  these  equations  are  solved  in  two  separate  phases 
so  that  fourth  order  characteristics  can  be  obtained  instead  of  the  cumber- 
some sixth  order  form.  In  the  first  stage  it  is  assumed  that  the  elevator  is 
fixed,  thereby  eliminating  the  elevator  equation  and  elevator  terms  from  the 
remaining  equations.  This  approach  affords  an  examination  of  the  short  period 
and  the  phugoid  modes.  The  second  stage  is  based  on  the  velocity  considered 
constant,  which  leads  to  the  exclusion  of  all  velocity  terms  and  a character- 
istic equation  representing  the  short  period  and  elevator  modes.  Neither  of 
these  simplifications  was  utilized  because,  particular  interest  was  centered 
on  the  phugoid  as  affected  by  the  free  elevator;  and  the  usually  accepted 
notion  of  the  phugoid  being  relatively  independent  of  the  elevator  motion  was 
not  believed  to  be  valid  for  systems  involving  fairly  strong  hinge  moment 
devices.  On  this  basis  a preliminary  analysis  was  undertaken,  applying  the 
more  complete  solution,  to  study  the  interdependence  of  the  various  derivatives 
and  enable  the  further  reduction  of  the  sixth  order  form.  A number  of  solu- 
tions were  made  for  the  modified  Navion  using  estimated  aerodynamic  character- 
istics and  considering  the  airplane  trimmed  in  level  flight  at  an  indicated 
airspeed  of  110  mph  and  an  altitude  of  8,000  feet*  Solutions  for  various 
static  stability  margins  were  obtained  both  analytically  and  by  an  analog 
computer.  In  the  computer  phase  of  the  work  the  analog  circuitry  duplicated 
the  revised  equations  as  given  and  combined  them  into  a single  dynamic  system. 
The  coefficients  of  the  characteristic  equation  required  for  the  analytical 
computations  were  determined  and  they  are  defined  as  follows: 
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A*  + AAV  CX3  y ZX  ^3=  c (28) 


/-?  — G>  ^ As  / 4 c,  A (^^We  '/U  3) 

3 = C^6(47t'^)^/4  ^ AA(C^-  C/>c,t//>3  C»WeA) 

aaT  ( c4/A  4 GX  G,^  ^ GG " at  (G&  * 4 O^g ) 

- G^  /G  1-  Ce. .4(  G~  GJ  (30) 

C = aa~  { GVG  ^4 1 Q ^ W)  *(/?,-  QA) 

" (4  +■  ~?'~  ( 9 ^ CyryA  Gg  ^ G O^c) 

- C^4  CV  y C^  (C^-iQi 

y G,  4/~ Q*  (Z^?c/&-  X G)4<G  ~ As  44 

- G C^Ue//A~  Cl  At  (Q  ' Gj  Q/s/4  ^ W) 

-C^/4  (31) 

z>  V /*  { G/ G (&+£}-(  4 y Q J7 

y G**  4/  y Co  Cajs)~  A G4)  ^ C 4 (G"  GJ) 
/C>G4  ^4)~  Os^  y G/6  AyJ^  A>cU  ( G G, 

* C^  C/4)  y C^(G6  y CG,r4  GV) 

y G**  G G6  y V6  { C 4 4a  QJ-  G * G,J] 

-Q<j}  (32) 

4~  =za  {4  ( G - G< ) { G^  /4  G0  _ G (4  y 0^)7 

* AUe  G c ] ~ A/z[CW8  (/?,  y G^,)^ 

~ G6  G^&7  ^ G^  ( C G-  ^ Qfc  C^)  ^ 

ZC  (4-Cu)-^  OJ-  C G^)  ^ A G^  GGe}  (33) 
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* c,  cO  c*.c)*> 


2.  Simplified  Solution 

An  examination  of  the  exact  solution  results  indicated  that  the  elevator 
mode  was  sufficiently  damped  to  permit  the  omission  of  the  damping  and  iner- 
tia terms,  and  /)  respectively,  from  the  hinge  moment  equation  with  a 

oC&  ' 

negligible  loss  in  accuracy.  This  condition  applied  even  when  bobweight 
systems  producing  large  hinge  moments  and  accompanying  increases  in  elevator 
inertia  were  considered.  With  the  elevator  equation  simplified  the  character- 
istic reduces  to  a fourth  order  which  can  be  solved  fairly  conveniently  by  a 
number  of  standard  mathematical  methods.  It  should  be  pointed  out  that  the 
simplification  eliminates  the  elevator  mode  but  does  include  the  effects  of 
the  free  elevator  on  the  remaining  short  period  and  phugoid  modes.  The  co- 
efficients of  the  quartic  in  terms  of  the  familiar  aerodynamic  quantities 
were  found  to  be 


A x t4-  /-) A3  /-  £ X1  C A ^ /)  — G 


(35) 


/?  - Q a C^/2  - — 


/ 


(36) 
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/ 


{<V  5?  ^ ' ‘Utfy  § ft' , c„J 

y-QlC^  +-  f^~(4  ~ QJ// 


(37) 


+ C^/t/cZ0-  ^U+6J]/-  CC, 

s)7  * 

c-  ( c4 


te/G 


7^ 


,.  _ d.*7£ 

A L_  «Ly  , — *■ 


*& 


} 


'"*cU  ~C~^  ( 

v*~  '/  £ 


(38) 


/O  = - 


l c:/r  _ u 

A - -S=*~  (4  ^4)  O 

_ iXs 


a. 


o 


X 


£ ~k  - £4,  G** 


c,  o. 


9 


(39) 


Effects  of  force  systems  are  accounted  for  by:  , the  hinge  moment  coeffi- 

c 

cient  due  to  devices  and  in  the  general  sense  that  resulting  as  well  from 
either  aerodynamic  causes  or  elevator  mass  unbalance  (effective  bobweight), 
/?'  v-:  the  elevator  mass  unbalance  as  acted  upon  by  the  airplane’s  linear 

acceleration,  and  $ , the  elevator  mass  imbalance  as  acted  upon  by  the  air- 
plane’s angular  acceleration.  Since  the  systems  under  study  are  essentially 
stick-free  static  stability  controls,  a rearrangement  of  the  coefficient 
equations  in  terms  of  familiar  static  parameters  and  simplified  dynamic  terms 
is  desirable.  This  would  facilitate  the  dynamic  analysis  by  clearly  indi- 
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eating  the  contributory  effects  of  the  static  stability  margins  and  by  re- 
ducing the  lengthy  mathematical  computations. 

Applying  the  relationships,  as  given  in  reference  1, 


C - - 

-4 

(40) 

4*^ 

r~ 

~'"Ue 

_ 

(4l) 

P 

r,  & 

(42) 

etc f 

(43) 

results  in  the  following  equality 


This  value  is  seen  to  be  a measure  of  the  aerodynamic  damping  of  the  airplane 


as  modified  by  the  free  elevator  and  is  referred  to  as  ^ ^ dx. 


is  identified  as  the 


£ 

In  a like  manner,  the  quantity  ^ 

effective  airplane  inertia  when  modified  by  a free  elevator  including  mass 
unbalance  and  is  denoted  by 

The  stick-free  level-f light  and  maneuver  margins  can  be  expressed  as  a 
function  of  repeated  groupings  of  related  terms  in  the  coefficient  equations. 
From  reference  1,  freeing  the  elevator  shifts  the  stick-fixed  neutral  point 
forward  as  given  by 


'X' -X*  = (M-Xce)- 


(45) 
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where 


(46) 


/l/  -X„_  = 


/V-  A/'- ^ 

yo  ' vo  /—  s’' 

C^< 


(47) 


Considering  the  dynamic  modes  at  essentially  unaccelerated  flight,  ( — , 

*•  t 

and  the  rearward  shift  in  neutral  point  caused  by  a bobweight  or  downspring 


can  be  developed  in  coefficient  form  as 

<4 


= - (a  a ) = C»'& 


-u  J 


/ 


UGi  G Cc  G, 

X c 


= c"-  C /- 


(48) 


(49) 


(50) 


c q5 

Thus,  the  general  expression  for  the  level-flight  stability  margin  may  be 
written  as  ^ 

* • ' - ' A C,  C>>;  L , x 

> -v-^-  (51) 

W . 


/y/-X  --  / ---  11^  y Gc 


C C 


At  ~Xcc,  = 


/ in  CGAf  . 

q,  ~cg1<52) 


The  rearward  displacement  of  the  stick-free  maneuver  point  from  the 
neutral  point  is  given  in  reference  1 as 

'= 


/k  - /v. 


-e  G . G g ) 


(53) 


are 


If  the  damping  contributed  by  components  other  than  the  horizontal  tail  is 

G~  \A//c 

neglected  and  the  relations  for  and  /v—  — 

<**  > W ^ c,  PC 

applied,  this  margin  becomes  ' y ^ 


(c. 


C,_ 


(54) 


The  shift-1  in  the  maneuver  point  due  to  bobweight  or  mass  unbalance  effects 
can  be  expressed  as 


/)( A/  ' ) =z  _ / 

1 

C,  7/' 

(55) 

following  equalities: 

(56) 

II 

ft 

(57) 

^ P ^ Ce 

(58) 

And  including  the  neutral  point  displacement  due  to  bobweights  or  downsprings 
as  previously  given,  the  following  generalized  representation  of  /CL  - /\/^  x 


is  obtained: 


^ ^ l Z Q ' 


(59) 


For  mass  balanced  elevator  systems  employing  a downspring  only,  ~ £)  and 
when  a bobweight  is  introduced  its  contribution  to  the  Ca  term  is  nullified 
by  the  /L/  term  as  shown.  Consequently  this  form  is  presented  to  retain  the 
general  meaning  of  CL  A as  the  sum  of  applied  hinge  moment  coefficients  from 
all  causes. 


The  stick-free  maneuver  margin  can  now  be  determined  from 

^(a'Za/o) 


(60) 


and  using  the  derived  relationships, 


/14  /hi?  - r q " 

n/c*-  fr  (4  +qjj> 


(61) 


By  substituting  all  the  foregoing  static  stability  and  damping  equiva- 
lences in  the  coefficient  equations  for  the  stability  quartic,  the  four  ex- 
pressions reduce  to 
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^ (Q  - Q<)  + c*&*.  + -jplc^  (M„-Xcf) 


(62) 


(63) 


C ( C~, 


)cC^ 


o^yj 


C /?'  l Cz>  O.^  (A/*,  /V rr7)  (li  Odx  (AL,  ~ A/  ) 


(64) 


^ ( ^'Uo  ^ 


D = 


/ 


y—2  x~ 
Lz  4_z . 


/f/>  Xcr) 


(65) 


x/  \ -2 

The  equations  can  be  somewhat  further  simplified  by  omitting  terms  which  are 
small  for  the  conditions  of  a particular  study.  The  selection  of  cruising 
or  high  speed  operation  for  an  airplane  with  an  unswept  moderate  aspect 
ratio  wing  will  permit  c to  be  neglected  relative  to  /X<L  ^ as  well  as 
the  product  of  and  other  comparitively  small  terms  (such  as  the  last 
quantity  in  the  equation  for^B).  In  addition,  as  a first  approximation  ^ 
may  be  considered  as  zero  for  balanced  elevators  with  bobweight  systems 


located  near  the  center  of  gravity  and  /)  may  be  neglected  relative  to 


A . 


These  changes  to  the  coefficient  equations  were  employed  in  the  subject 
investigation.  Power  effects  may  be  taken  into  account  by  including  the 


c 


*7 


Zu 


term  in  the  pitching  moment  equation  of  motion.  This  results  in  the 


addition  of  the  increment  A ( A/c  y 


c. 


1U 


-2.  Cc 


to  the  expression  for  the 


level-flight  stability  margin.  As  the  maneuvering  stability  applies  to  con- 
stant velocity  operation  it  is  unaffected  by  power  variations. 


3.  Theoretical  Results 

An  extensive  computational  program  was  completed  using  the  simplified 
solution  to  obtain  the  short  period  and  phugoid  mode  characteristics  as  a 
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function  of  the  independent  variables, /I X c<$  and  — SVC‘  . With 

these  parameters  as  the  prime  criteria  for  the  theoretical  analysis  the  ef- 
fects of  center  of  gravity  displacement,  bobweights  and  downsprings  can  be 
isolated  for  study  in  a manner  consistent  with  the  objectives  of  the  test 
program.  As  has  been  demonstrated  for  a given  airplane,  the  maneuver  margin 
can  be  controlled  by  a combination  of  center  of  gravity  location  and  bob- 
weight  strength  and  the  difference  between  the  maneuver  and  neutral  points 
can  be  controlled  solely  by  a downspring.  For  convenience,  the  increment 

- /Vc'  is  called  the  static  damping  parameter  since  this  quantity  is 
dependent  upon  the  aerodynamic  damping  in  pitch  for  airplanes  without  control 
devices  and  it  also  is  an  important  factor  affecting  the  damping  of  the  dy- 
namic characteristics. 

A single  operating  condition  was  considered  for  the  modified  Navion  air- 
plane. This  was  trimmed  level  flight  at  an  indicated  airspeed  of  110  mph 
and  an  altitude  of  8,000  feet  with  the  gross  weight  assumed  constant  at  2,500 
lbs.  Values  of  /)/■>„  and  A/^  ~ /V0 ' were  selected  (thus  fixing 

AV o ~ A(  c$  ) covering  a wide  range  of  stability  margins  and  the  coeffi- 
cients of  the  characteristic  equation  were  then  determined  using  estimated 
values  of  the  other  aerodynamic  terms  in  the  simplified  equations.  The  ap- 
proximations for  the  aerodynamic  and  physical  parameters  applied  in  this 
analysis  are  given  in  Appendix  II.  Solutions  for  the  roots  of  approximately 
eighty  stability  quartics  involved  were  obtained  by  a semi-empirical  method 
which  allowed  a good  degree  of  accuracy  without  incurring  any  unreasonably 
lengthy  computational  work.  Conversions  were  made  of  negative  real  roots  and 
negative  real  parts  of  complex  roots  to  the  reciprocal  of  the  time  to  damp 
to  half  amplitude,  positive  real  roots  and  positive  real  parts  to  the  recip- 
rocal of  the  time  to  diverge  to  double  amplitude,  and  the  imaginary  parts  of 
complex  roots  to  oscillatory  period.  Inverse  quantities  are  used,  similar  to 
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root  loci,  so  that  the  damping  characteristics  can  he  plotted  as  continuous 
functions  through  the  neutral  stability  boundary,  thereby  simplifying  the 
definition  of  all  trends.  The  complete  theoretical  results  are  presented  in 
Figure  4 in  the  form  of  damping  and  period  as  a function  of  c<$  for 

various  fixed  values  of  /V^  ~~  A/0  . The  stick-free  maneuver  margin  range 

extends  from  large  positive  or  highly  stable  values  to  large  negative  or 
highly  unstable  ones.  The  magnitude  of  the  static  damping  parameter  varies 
from  .04,  comparable  to  the  modified  Navion  configuration  without  downspring, 
to  -.20,  that  obtained  with  a very  strong  downspring.  Oscillatory  motion, 
determined  by  real  parts  of  complex  roots,  is  identified  by  dotted  lines  and 
aperiodic  motion,  determined  from  real  roots,  is  identified  by  solid  lines. 

An  examination  of  the  damping  characteristics  in  the  upper  part  of  the 
plot  reveals  a number  of  significant  facts.  First,  the  short  period  mode 
has  no  importance  in  the  subject  study.  In  almost  the  entire  stable  maneuver 
margin  range  this  mode  is  independent  of  both  parameters  and  maintains  a 
constant  damping  characteristic  of  the  Navion.  At  maneuver  margins  less  than 
Xc$  = • 048 , it  is  shown  as  a convergent  aperiodic  but  actually  has 
no  bearing  on  the  resultant  airplane  motion  because  the  phugoid  mode  is  then 
moderately  or  highly  divergent  and  thus  it  has  the  predominating  effect. 
Second,  the  suspicions  that  prompted  the  program  are  substantiated  since  both 
bobweights  and  downsprings  have  a considerable  influence  on  airplane  dynamics 
This  is  made  evident  by  the  phugoid  curve  family  which  rapidly  deteriorates 
from  the  normal  damped  oscillation  at  moderately  stable  maneuver  margins  ’to 
highly  unstable  oscillations  and  eventually  pure  divergences  at  reduced 
levels  of  stability.  And  finally,  the  quantity  //^r  A/0  is  seen  to  represent 
a generalized  parameter  of  importance  in  the  study  of  dynamic  stability 
effects.  Of  only  academic  interest  are  the  indications  that  the  short  period 
and  phugoid  modes  are  joined  by  two  aperiodic  branches  which  subsequently 
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develop  into  a single  oscillatory  mode.  Tnese  motions  are  completely  domi- 
nated, as  is  the  convergent  short  period  branch,  by  the  concurrent  highly  di- 
vergent phugoid. 

The  period  trends  shown  in  the  lower  part  of  the  plot  correspond  as 
would  be  expected  to  the  associated  damping  curves  --  in  the  sense  that  as 
the  oscillation  becomes  increasingly  divergent  the  period  increases,  asymp- 
toticaly  approaching  infinity  at  the  point  where  the  aperiodic  motion  origi- 
nates. With  reference  to  the  normal  short  period  mode,  decreasing  the  man- 
euver margin  increases  the  period  and  variations  in  the  static  damping 
parameter  alter  the  phugoid  slightly. 

It  is  evident  from  the  results  presented  in  Figure  4 that  the  correction 
of  rearward  center  of  gravity  locations  by  progressively  applying  bobweights 
has  a continuously  favorable  effect  on  the  phugoid  damping.  Whereas  when 
downsprings  are  introduced  the  damping  is  improved  until  some  optimum  spring 
strength  is  reached  and  then  additional  increments  of  downspring  cause  the 
damping  to  become  progressively  reduced.  The  overall  trends  indicate  that 
the  bobweight  is  a considerably  more  powerful  device  in  the  control  of  the 
dynamic  characteristics  than  is  the  downspring.  This  is  clearly  borne  out  by 
the  fact  that  as  the  maneuver  margin  is  made  less  stable  the  phugoid  rapidly 
increases  in  divergence.  Changes  in  the  static  damping  parameter  displace 
these  variations  but  do  not  alter  the  general  nature  of  the  unstable  trend. 

In  order  to  more  conveniently  isolate  these  effects  with  reference  to  stabil- 
ity boundaries  and  correlate  them  with  the  flight  data,  crossplots  of  the 
phugoid  curves  were  made  with  lines  of  constant  damping  times  and  period  con- 
structed  on  a /Vyy,~  A/0  versus  /4£,  ~ Xc$  grid.  The  data  were  obtained  by 
interpolation  and  are  limited  to  the  portion  of  the  complete  plot  which  could 
practically  be  explored  in  the  flight  program.  The  damping  crossplot  is 
given  in  Figure  5 and  the  period  crossplot  is  given  in  Figure  6. 
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Three  separate,  fairly  well  defined  stability  regions  were  established 
as  seen  in  Figure  5*  These  consist  of  aperiodic  divergence,  damped  oscilla- 
tions and  undamped  or  divergent  oscillations.  The  important  fact  indicated 
by  this  plot  is  that  if  a damped  or  essentially  neutrally  stable  phugoid  is 
to  be  maintained  the  maneuver  margin  must  be  stable  and  the  static  damping 
parameter  must  be  positive  (unless  the  maneuvering  stability  is  moderately 
high  in  which  case  sufficient  spring  may  be  introduced  to  cause  the  damping 
parameter  to  become  negative).  It  can  be  seen  that  increasing  the  bobweight 
strength  or  positive  /\ ^ c$  values  displaces  the  phugoid  from  the  pure 
divergent  region  to  either  the  damped  oscillatory  one,  at  positive  /%? 
values  representing  little  or  no  downspring,  or  to  the  undamped  oscillatory 
one,  at  negative  /c  values  representing  moderate  or  heavy  downsprings. 

Increasing  the  downspring  magnitude  when  the  maneuver  margin  is  moderately 
or  highly  stable  displaces  the  phugoid  from  a damped  oscillatory  region  to 
an  undamped  one.  At  low  positive  maneuver  margins  adding  downspring  dis- 
places the  phugoid  from  the  divergent  aperiodic  region  to  first  the  damped 
oscillatory  one  and  finally  to  the  undamped  oscillatory  region.  At  negative 
maneuver  margins  the  displacement  is  from  the  divergent  aperiodic  region  to 
the  undamped  oscillatory  region. 

The  period  crossplot  of  Figure  6 shows  that  the  trends  of  the  constant 
period  lines  are  consistent  with  that  of  the  aperiodic  boundary  as  would  be 
expected.  In  the  oscillatory  region  considered  the  period  varies  from  12  to 
60  seconds. 

4.  Analog  Computer  Studies 

A number  of  theoretical  solutions  of  the  modified  Navion  dynamics  were 
also  obtained  by  analog  computer  in  order  to  observe  the  combined-mode 
transient  response  of  the  airplane,  including  any  peculiarities  of  the  over- 
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all  motion  not  evident  in  the  analytical  treatment,  as  well  as  to  make  spot 
checks  of  the  responses  measured  in  the  test  phase.  For  this  purpose  aero- 
dynamic derivatives  extracted  from  static  flight  data  were  used  in  conjunction 
with  the  complete  equations  of  motion.  The  experimentally  derived  quantities 
and  the  reduction  methods  applied  are  presented  in  Appendix  III.  Phugoid 
responses  to  constant  elevator  inputs  were  obtained  for  various  force  device 
arrangements  and  the  same  operating  conditions  were  considered  as  in  the 
analytical  analysis.  The  results  for  six  different  dynamic  runs,  transformed 
into  comparable  inputs,  are  presented  in  the  summary  of  flight  data  for  com- 
parative study. 
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VIII.  FLIGHT  TEST  PROCEDURE 


A.  Static  Testing 

A series  of  static  stability  tests  were  conducted  at  various  center  of 
gravity  positions  with  different  magnitudes  of  bobweight  and  downspring  in 
order  to  determine  the  stick-free  static  stability  characteristics  required 
for  the  dynamic  analysis.  Force  gradients  were  measured  for  both  level-flight 
and  maneuvering  flight  to  establish  the  effectiveness  of  the  force  devices  in 
the  control  of  the  stick-free  neutral  and  maneuver  point  locations  For  the 
sake  of  completeness  elevator  deflection  variations  were  also  measured  to 
obtain  the s tick-fixed  neutral  and  maneuver  point  locations.  All  straight 
and  level  runs  were  made  at  indicated  airspeeds  ranging  from  80  to  ll|0  mph 
with  the  control  force  trimmed  at  110  mph.  The  maneuvering  tests  were  made 
in  steady  accelerated  turns  flown  at  normal  load  factors  ranging  from  1.2 
"g"s  to  2.0  "g"s  with  the  indicated  airspeed  held  constant  at  110  mph  and  the 
control  force  trimmed  for  straight  and  level  flight.  An  average  pressure 
altitude  of  5>000  feet  was  maintained  during  all  testing. 

Speed  points  and  accelerated  turns  were  carefully  stabilized,  gradients 
were  given  rough  checks  in  flight  and  repeat  runs  were  made  where  necessary. 
Air, speed,  normal  load  factor  and  indicator  readings  of  elevator  deflection 
and  stick  force  were  visually  recorded  during  the  tests. 

Preliminary  flight  checks  were  made  to  establish  the  constancy  of  the 
bobweight  and  downspring  hinge  moments  over  the  test  range  of  stick  deflec- 
tion. It  was  also  confirmed  that  the  downspring  did  not  affect  the  maneuver- 
ing stability.  Sufficient  data  were  then  taken  to  accurately  locate:  the 

neutral  and  maneuver  points  for  the  balanced  elevator  without  devices,  the 
neutral  points  for  bobweight  and  downspring  combinations  inducing  up  to  32 
lbs.  of  equivalent  stick  force,  and  the  maneuver  points  for  bobweight 
strengths  up  to  2k  lbs.,  all  in  control  force  increments  of  8 lbs.  each. 
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B.  Dynamic  Testing 


The  general  method  of  dynamic  testing  employed  is  that  duplicating  as 
closely  as  possible  the  theoretical  case  of  transient  response  to  a step  in- 
put in  elevator  control.  This  was  done  by  trimming  the  airplane  at  an  indi- 
cated airspeed  of  110  mph  in  straight  and  level  flight  at  5>000  ft.,  then 
applying  a sufficient  steady  control  force  to  attain  equilibrium  at  a speed 
between  2 and  5 mph  above  or  below  the  trim  velocity  and  then  instantaneously 
releasing  the  stick.  Airspeed  time  histories  were  obtained  by  photorecorder 
for  the  full  duration  of  each  phugoid  response.  Considerable  care  was  taken 
to  achieve  the  trim  speed  accurately  and  avoid  perceptible  pitch  inputs  and 
gust  disturbances  so  that  true  responses  could  be  recorded. 

In  order  to  substantiate  the  significant  theoretical  trends,  it  was 
planned  to  explore  as  great  a portion  of  the  static  stability'  grid  of 
Figure  5 as  possible.  The  actual  test  conditions  were  somewhat  limited 
by  several  practical  considerations.  Center  of  gravity  locations  were  re- 
stricted to  positions  forward  of  36.3%  MC  due  to  a critical  shortage  of 
down  elevator  deflection  required  for  recovery  frcm  rapid  divergences.  The 
positions  for  the  dyrnamic  tests  were  at  approximately  32%,  3l$  and  36%  MAC. 

In  addition,  the  maximum  effectiveness  of  the  modified  elevator  tab  limited 
the  total  control  force  imposed  by  the  devices  to  32  lbs.  Individual  bob- 
weight  and  downspring  forces  were  varied  from  0 to  2h  lbs.  and  numerous 
phugoid  responses  were  recorded  for  various  combinations  of  each  at  the 
three  center  of  gravity  locations. 
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IX.  FLIGHT  DATA  ANALYSIS 


A.  Static  Stability 

Indicator  readings  of  stick  force  and  elevator  deflection  were  trans- 
formed to  actual  values  using  current  calibrations.  Lift  coefficients  and 
free  stream  dynamic  pressures  were  computed  for  the  level-flight  runs  from 
the  take-off  weight  and  the  indicated  airspeed  readings.  Plots  were  then  con- 
structed of:  versus  CL L for  each  combination  of  center  of  gravity 

and  total  control  force  due  to  the  bobweight  and  downspring,  versus 

??  for  each  combination  of  center  of  gravity  and  bobweight  force,  and  (5 
versus  (L L & versus  /V  for  each  center  of  gravity.  The  slopes  of  all 

the  faired  variations  or  gradients  representing  a measure  of  the  static 
stability,  were  then  measured.  Typical  flight  data  are  shown  in  Figure  7. 

It  can  be  seen  that  linearity  is  obtained  in  the  stick-fixed  curves  and  the 
maneuver  force  gradient.  With  regard  to  the  level-flight  force  character- 
istics some  variation  in  C2a  with  elevator  deflection  is  evident.  The 
deviation  from  the  straight  slope  for  other  tests  varied  more  or  less  from 
the  example  presented  depending  upon  the  particular  range  of  elevator  deflec- 
tions involved.  By  taking  the  slope  of  the  /~$  /(^  curve  at  trim  this  tend- 
ency does  not  impair  the  accuracy  of  the  stability  measurement. 

All  force  and  elevator  gradients  were  compiled  in  summary  plots  as  a 
function  of  center  of  gravity  location.  These  results  are  given  in  Figure 
8.  The  stick-fixed  and  stick-free  neutral  and  maneuver  points  are  indicated 
at  the  zero  gradient  intercepts.  To  aid  in  the  analysis  the  curves  were 
faired  in  accordance  with  the  required  relationship  between  the  slopes  of 

and  c£ /c/r) versus  center  of  gravity.  A similar 
procedure  was  followed  as  a guide  in  fairing  the  test  points  for  the  case  of 
atS/rAC,  and  as  a function  of  center  of  gravity.  It  can  be 

shown  that 
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where  -^-  , the  average  wing  loading  is  15.2  lbs. /ft. ^ and  , the  trimmed 

lift  coefficient,  is  0.^9.  The  overall  correlation  of  the  static  results  is 
satisfactory;  showing  relatively  little  data  scatter,  linear  stability  grad- 
ients and  agreement  between  the  neutral  and  maneuver  point  displacements 
caused  by  equal  increments  of  device  force.  The  static  stability  character- 
istics are  therefore  considered  as  having  been  accurately  determined. 


B.  Dynamic  Stability 

Approximately  eighty  phugoid  runs  were  recorded  for  twenty-nine  different 
combinations  of  maneuver  margin  and  static  damping  parameter.  Faired  enve- 
lopes of  each  oscillatory  response  were  constructed  in  order  to  qualitatively 
evaluate  the  damping  and  period  characteristics.  It  was  evident  that  gust 
inputs  during  the  course  of  the  motion  had  caused  discrepancies  among  the 
various  records  for  the  same  test  conditions.  A judicious  selection  of  the 
most  reasonable  indications  was  made  on  the  basis  of  a cyclic  consistency 
approaching  a constant  period  and  a continuous  change  in  amplitude.  Aperiodic 
motions  being  of  short  duration  were  not  as  subject  to  gust  effects  and 
showed  satisfactory  repeatability.  The  velocity  time  histories  selected  for 
each  condition  of  static  stability  are  presented  in  Figures  9 through  15. 

For  comparison  purposes  most  of  the  original  data  were  proportionately  trans- 
formed into  velocity  variations  for  a positive  input  of  5 mph.  The  actual 
mean  velocity  differed  slightly  from  the  110  mph  value  due  to  the  residual 
friction  level  and  the  difficulty  in  trimming  when  the  force  gradients  are 
low . 
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Dynamic  results  have  teen  arranged  in  static  stability  groupings  to 
facilitate  the  correlation  of  the  flight  data  with  the  theoretical  trends 
previously  described.  The  dynamic  equivalence  of  bobweight.  and  center  of 
gravity  effects  on  the  maneuver  margin  is  first  demonstrated  by  a series  of 
phugoid  modes  shown  in  Figure  9*  In  this  study  different  combinations  of 
bobweight  and  center  of  gravity  giving  approximately  the  same  maneuver  margin 
are  compared  for  various  static  damping  parameters.  It  can  be  seen  that  for 
a given  value  of  downspring,  equivalent  modes  are  obtained  when  the  maneuver 
margins  are  comparable  regardless  of  the  bobweight  magnitude  and  center  of 
gravity  location.  The  discrepancies  shown  are  as  expected  considering  the 
actual  differences  in  maneuver  margins  in  each  pair.  In  Figures  10  through 
12  the  static  damping  parameter  is  held  constant  at  various  values  and  the 
maneuver  margins  are  progressively  reduced.  It  is  evident  that  reductions 
in  the  maneuvering  stability  have  continuously  deleterious  effects  on  the 
phugoid  mode.  This  trend  is  unaffected  by  the  static  damping  parameter  as 
shown  for  the  zero  downspring  force  in  Figure  10,  the  moderate  spring  in 
Figure  11  and  the  heavy  spring  in  Figure  12.  In  Figures  13  through  15  the 
maneuver  margins  are  maintained  approximately  constant  at  several  levels  and 
the  static  damping  parameter  is  progressively  reduced  by  steadily  increasing 
the  downspring  force.  In  Figure  13  it  is  indicated  that  when  the  maneuver 
margin  is  high  the  introduction  of  downspring  forces  has  an  adverse  effect, 
changing  the  phugoid  from  a damped  oscillation  to  an  undamped  oscillation. 
Figures  Ik  and  15  demonstrate  the  improvement  in  the  dynamics  attainable  at 
low  positive  maneuver  margins  with  a moderate  application  of  downspring. 

It  is  apparent  from  the  grouping  in  Figure  14  that  further  increases  in 
downspring  force  have  the  effect  of  partially  nullifying  the  original  gains. 
On  the  basis  of  the  foregoing  results  the  theoretical  effects  of  force  de- 
vices on  the  dynamics  are  considered  as  having  been  qualitatively  confirmed 
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in  the  flight  test  phase. 

For  general  interest  the  damping  time  and  period  of  each  phugoid  response 
was  measured  and  a quantitative  comparison  with  the  theoretical  results  was 
made.  These  results  are  included  in  Figures  5 and  6.  This  presentation 
offers  an  alternate  indication  of  the  qualitative  agreement  with  theory  and 
also  shows  an  unexpected  agreement  in  a quantitative  sense.  It  was  originally 
thought  that  the  necessary  theoretical  assumption  of  motion  limited  to  small 
perturbations  and  purely  linear  aerodynamic  effects  would  make  a quantitative 
check  unlikely.  The  analog  computer  comparisons  that  are  included  in  Figures 
10,  11,  12  and  Ik  are  also  satisfactory  considering  these  same  factors. 
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X*  CONCLUSIONS 


When  auxiliary  control  force  devices  such  as  bobweights  and  downsprings 
are  used  to  obtain  satisfactory  static  stability,  consideration  should  be 
given  to  the  resultant  effects  on  the  phugoid  mode*  Conversely,  it  may  be 
possible  to  correct  inherently  undesirable  phugoid  characteristics  by  a 
judicious  application  of  force  devices* 

A bobweight  will  have  no  adverse  influence  on  the  phugoid  and  is  likely 
to  produce  appreciable  improvement  in  the  damping*  A limited  use  of  down- 
spring may  have  favorable  effects  on  the  phugoid  but  an  excessive  application 
will  induce  strong  instabilities.  On  the  dynamic  basis  the  bobweight  is 
preferable  to  the  downspring  as  a static  stability  control  device. 

After  considerable  flight  operation  with  the  center  of  gravity  located 
aft  of  the  stick-fixed  neutral  point  the  pilot  concluded  that  no  adverse  ef- 
fects due  to  the  static  instability  were  noticeable. 
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XI.  RECOMMENDATIONS 


To  enlarge  the  scope  of  the  general  study  of  the  phugoid  it  is  recom- 
mended that  stabilization  means  such  as  pitch  damping  mechanisms  be  investi- 
gated. It  would  also  be  advisable  to  extend  the  present  study  to  include 
flight  conditions  other  than  cruising  operation.  In  order  to  increase  the 
practicality  of  the  results  derived  from  phugoid  studies  it  is  recommended 
that  minimum  requirements  for  satisfactory  flying  qualities  be  established 
for  this  mode. 
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XIII.  NOMENCLATURE 
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Ce 

C* 
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Q 

C 

G, 
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4 
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& 

/? 
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4 

///? 
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4 

4 

?/? 

We 

4% 
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mean  aerodynamic  chord  of  wing,  ft. 
root -.mean- square  chord  of  elevator,  ft. 
root-mean- square  chord  of  elevator  trim  tab,  ft. 
airplane  drag  coefficient 
elevator  hinge  moment  coefficient 
airplane  lift  coefficient 

pitching  moment  coefficient  about  airplane  center  of  gravity 

applied  hinge  moment  coefficient 

stick  control  force,  positive  for  push,  lbs. 

stick  force  equivalent  to  bobweight  hinge  moment,  lbs. 

stick  force  equivalent  to  downspring  hinge  moment,  lbs. 

acceleration  due  to  gravity  (32.2  ft. /sec  .) 

airplane  inertia  term 


elevator  mass  unbalance  effect 


elevator  inertia  term  — ) 

elevator  hinge  moment,  positive  when  deflecting  the  trailing  edge 
downward , lb . - f t . 

elevator  radius  of  gyration  about  the  hinge  line,  ft. 
airplane  radius  of  gyration  about  the  pitch  axis,  ft. 
elevator  mass  unbalance  term  as  affected  by  airplane  pitching  ac- 
celeration ( j 

( ^er*c*  J 


horizontal  tail  moment  arm,  ft. 
mass  of  airplane,  slugs 
mass  of  elevator,  slugs 

stick-fixed  neutral  point  location  in  terms  of  (S 
stick-free  neutral  point  location  in  terms  of 
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/14 

/k; 

77 

9 

•vSe 

£ 

2d 

1/ 

Alt 

u/ 

% 

<=K_ 

^-t: 

6 

6^ 

£ 

© 
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stick-fixed  maneuver  point  location  in  terms  of  CZ 

stick-free  maneuver  point  location  in  terms  of  (2 

normal  load  factor,  acceleration  in  terms  of  ^ 

period  of  oscillation,  sec. 

free  stream  dynamic  pressure,  lbs. /ft. 

2 

wing  area,  ft. 

2 

elevator  area,  ft. 

2 

elevator  tab  area,  ft. 

time  required  for  oscillation  to  damp  to  half  amplitude,  sec. 
time  required  for  oscillation  to  double  amplitude,  sec. 
time,  sec. 

flight  velocity  ratio 
true  flight  velocity,  ft. /sec. 

change  in  true  velocity  from  trimmed  value,  ft. /sec. 
calibrated  airspeed,  ft. /sec. 

change  in  calibrated  airspeed  from  trimmed  value,  ft. /sec. 

airplane  gross  weight,  lbs. 

2 

wing  loading,  lb 3. /ft. 

center  of  gravity  location  in  terms  of  CZ 

distance  from  hinge  line  to  center  of  gravity  of  equivalent  elevator 
mass,  positive  when  aft  of  hinge  line,  ft. 

airframe  angle  of  attack,  rad. 

horizontal  tail  angle  of  attack,  rad. 

elevator  deflection,  positive  for  trailing  edge  down,  rad. 
trim  tab  deflection,  positive  for  trailing  edge  down,  rad. 
angle  of  downwash  over  horizontal  tail,  rad. 
angle  of  pitch  of  airplane,  rad. 
root  of  stability  equation 
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yU-o, 

? 

T: 


relative  density  term  of  airplane 
relative  density  term  of  elevator 


free  stream  air  density,  slugs/ft. 


f -3C 

/ 'Ve- 

( 


time  factor 


(V^V) 


elevator  effectiveness 


Whenever  2^ 


^ . 5 


ct  6 

a k 


subscripts,  a derivative  is  indicated.  For*  example,  ^ . 

C =-^ 


and £<4 are  us t 

<3  °<L 


c)  ) 

y]/^- /y/  static  damping  parameter. 


d as 
and 
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APPENDIX  I 

PHYSICAL  CHARACTERISTICS  OF  MODIFIED  NAVION  AIRPLANE 

Wing 


Area 

184.2  ft. 

Mean  Aerodynamic  Chord 

5.70  ft. 

Aspect  Ratio 

6.04 

Span 

33.36  ft. 

Root  Chord 

7.22  ft. 

Tip  Chord 

3.83  ft. 

Airfoil  Section 

Root 

NACA  4415 

Tip 

NACA  6410 

Horizontal  Tail 

Area 

35.2  ft. 

Mean  Chord 

3.48  ft. 

Aspect  Ratio 

2.96 

Span 

10.22  ft. 

Root  Chord 

4.00  ft. 

Tip  Chord 

2.89  ft. 

Tail  Moment  Arm 

14.93  ft. 

Elevator 

Area 

12.3  ft.‘ 

Root -Mean -Square  Chord 

1.23  ft. 

Span 

10.22  ft. 

Root  Chord 

1.40  ft. 

Tip  Chord 

1.03  ft. 

Aerodynamic  Balance 

none 
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APPENDIX  II 


ESTIMATED  PARAMETERS  APPLIED  IN  THE  ANALYTICAL  ANALYSIS 


Operating  Conditions 


)/c  = 110  mph 
Alt.  = 8,000  ft. 
IV  = 2,500  Lbs. 


C = *556 


= .0443 

C = 4.58 

rad.  1 

G*  = .326 

rad.  ^ 

= -.674 

rad.  1 

(Z,^  = -.0818 

rad. 

-.0312 

rad.  1 

Cj,  = -.206 

rad.  1 

C/;h  = -.607 

rad.  ^ 

^ = .55 

.025 

0 

0 


1+1 
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APPENDIX  III 


EXPERIMENTAL  DERIVATION  OF  PARAMETERS  APPLIED  IN  ANALOG  COMPUTER  STUDIES 


Operating  Conditions 


V = 110  rrrph 
Alt.  = 5,000  ft. 
]V  = 2,800  lbs. 


■.  O- 


Based  on  the  wing  planform  and  section  characteristics  the  airplane  lift 


curve  slope  is  estimated  as 


(t  = 4.58  rad. 


The  trimmed  lift  coefficient  is 


f 

^c>  = 15.2  lbs. /ft. 


= (no) 

" 391 

= 31.0  lbs. /ft.' 


(III-1) 


r 15^2 

v*  _ 31.0 


= ,49 


C_.C 


The  drag  coefficient  derivative  is  given  by 


Q,  = Q,  * c. 


(III-2) 


The  airplane  drag  equation  as  determined  from  gliding  flight  tests  at  wind- 
milling  power  and  from  the  estimated  propeller  characteristics  is 


= .0290  + .0605  Ct 

Q = .121  CL 


(m-3) 


(m-4) 


WADC  TR  53-350 


CT/>  = .121  x .i| $ x 1+.58 
=*< 

= .272  rad."1 

The  drag  coefficient  considered  in  the  equations  of  motion  represents 
the  variation  of  force  along  the  flight  axis  with  velocity  and  consequently 
the  contribution  of  thrust  changes  with  airspeed  should  be  included. 

Effective  Q = Q - ^gLZl  „ —J—  (IH-5) 

Q z -0290  i-  .0605  (. h9 )2 
- .0435 

Power-On  glide  and  climb  tests  were  conducted  to  establish  the  thrust 
derivative  term  for  the  given  operating  conditions. 


/ 

^ J/  |/ 

Effective  C2 . 


-.0120 

.01+35  /-  .0120 
.0555 


The  pitching  moment  derivative  'with  respect  to  angle  of  attack  is 
computed  from  the  stick-fixed  neutral  point  for  a specified  center  of 
gravity  location. 

(A i -Xce)  (III-6) 

= -I+.58  (.308  - Xc$) 

- 4.58  -l.l-lio  (rad."-*-) 


The  elevator  effectiveness  is  determined  from  the  stick-fixed  static 
stability  tests. 

Level  Flight  Data 


o - - 57.3 

^8  ~ dXcQ 


(III-7) 
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1+3 


o£*& 

cXCi.  dXcL 


58.6  deg. 


c - . Shi 

tr«  £ 


Maneuvering  Flight  Data 


c = 


58.6 


- .978  rad.' 


37.3 


c£*6 


c/7?  ctXco 


'X  "}?  cX.  Xec? 

29.3  deg. 


(III-8) 


C _ . 57.3  x ,49 

^c"  2Q.1 


29.3 

= -.958  rad. 


-1 


Average 


G, 


G*  = ~ *968  rad. 


-1 


The  elevator  hinge  moment  variation  with  elevator  deflection  is  deter- 
mined from  the  stick-free  static  stability  tests. 


Level  Flight  Data 


G 


Cw±  „ d^C/3i 


Q -5e  Ce  yt  ^ Ct  <XXc$ 

C?  ) the  control  system  gearing  ratio  = 1.05  ft.  ^ 
^ the  tail  efficiency  =0.9 

_5e=  12.3  ft.2 
C-  = 1-23  ft. 


(111-9) 


= 20.2  ft.2 

.968  x 20.2 

1.05  x 12.3  x 1.23  x 0.9 
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IA- 


Maneuvering  Flight  Data 


(C  * - 1.370  rad. 


G = 


s ctr)c£Xc9 


(111-10) 


c£jf± 

c£  >7  o£Xc$ 


= 307.5  lbs. 


Q- 


- .968  x 307.5 

15.2  x 1.05  x 12.3  x 1.23  x 0.9 


= -I.370  rad. 


The  elevator  hinge  moment  variation  with  angle  of  attack  was  obtained 
from  the  measured  difference  between  the  stick-fixed  and  the  stick-free  neutral 


points. 


Cj  = - — (a/0 -a/A  (in-11) 

^ \ ' 


Ai~A/0  = .308  -.301 
= .007 


G = - 


4.58  x 1.370  x .007 

= " .968 

= -.0454  rad.  1 


The  damping  in  pitch  was  determined  by  measuring  the  difference  in  ele- 
vator deflection  required  for  trim  between  level  flight  operation  and  steady 
maneuvering  flight  at  99  =1*5*  Incremental  elevator  deflections  at  various 
lift  coefficients  were  then  converted  to  pitching  moment  increments  at  given 
pitch  velocities  by  applying 


<4  CL  - CL.  * a <5 


(m-12) 
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(III-13) 


# 

where  = the  pitching  velocity,  rad. /sec. 

^ = .00205  slugs/ft.3 

A plot  Of  AC- yy,  as  a function  of  GL  was  constructed  yielding  a linear 

tion  with  a slope  of  = -.119  sec. /rad.  And  in  non-dimensional  form  the 

clfr 

damping  derivative  is 


vana- 


r - 
v-^  - 


^2. 


/ 


a&  ot  IT 

= 1.32  sec. 


(III-14) 


c 

cle 


ai2 

1.32 


= -.0902  rad 


-1 


(~>rWc* 


The  damping  due  to  downwash  lag  is  given  hy 


C -C 

^—>7?  . Vs — >T) 


cA_  <5 


(m-15) 


The  downwash  derivative  is  estimated  to  he 

- = .382 

aot 

= -.0902  x .382 

C&< 

= -.0345  rad.-1 


The  change  in  hinge  moment  with  pitching  velocity  is  computed  from 


(UI-16) 

G = 

/ — cA(f 

(III-17) 

^ = 14.93  ft. 
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¥> 


yU  = 40.9 

C = 5.70  ft. 


Ca  . = 


.0454  x 14.93 
I-.382)  40.9x5.70 


= -.00470  rad. 


The  hinge  moment  due  to  downwash  lag  is  obtained  from 


Q = Q 

cAc>{_  de 

m -.00470  x .382 


(111-18) 


= -.00180  rad. 


The  elevator  damping  derivative  is  calculated  theoretically  by  the 
method  presented  in  reference  1. 

C^s=  [C  *■  £>(Cc^)±  ] (m-19) 

Q , the  mean  chord  of  the  horizontal  tail  = 3.48  ft. 

t the  estimated  horizontal  tail  lift  curve  slope  = 2.58  rad.  ^ 
and £) , constants  dependent  upon  the  chordwise  hinge  line  location,  are 
1.35  and  .08,  respectively. 


'^cte  ~ 2x4o 


i^TO  6-35  + -08  X 2.58) 


= - .0116  rad. 


The  applied  hinge  moment  coefficient  due  to  force  devices  is  given  by 


c = 


b G ft  Ce 



"1,05  x 31.0  x 12.3  x 1.23 

■ •oozos 


(III-20) 
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The  airplane  inertia  parameter  was  obtained  from  the  measured  moment  of 
inertia  about  the  jack  points.  After  applying  corrections  for  the  actual 
test  loading  conditions  and  the  transfer  to  the  center  of  gravity,  the  moment 
of  inertia  was  determined  as  = 3,800  slug-ft.2 


4 = ^7# 


-Zc<£ 


(III-21) 


2 x 32.-2  x 3.800 
= 2,800  x 40.9  (5.70)** 

= .O658 


The  effective  elevator  inertia  parameter  was  determined  from  estimations 
of  the  basic  elevator  moment  of  inertia  and  the  contribution  of  the  bob- 
weight  installation  including  the  constant  increment  required  to  initially 
balance  the  control  system. 

2 

Basic  I = 0.2  slug-ft. 

Neg.  Bobweight  I = 0.8  slug-ft.2 

Pos.  Bobweight  I = 0.2  (slug-ft.2) 


= 0.2  + 0.8  + 0.2  /£#  (III-22) 

-.1.0  + 0.2  /J*,  (slug-ft.2) 

^ = J 3 Tctycf-  (III-23) 

2 (1.0  + 0.2  ) 

.00205  x 12.3  x 1.23  (40.9  x 5.70)^ 

= .0012  + .00024 


The  elevator  mass  unbalance  term  is  computed  from  the  bobweight  effect. 

A = — (iii-24) 

' ^88  "“*e  C-GyA6  c 

h B 
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Figure  1.  Modified  Navion  Airplane,  General  View  of  Test  Airplane 


Original  Tail 


Figure  2.  Modified  Navion  Airplane,  Modified  Horizontal  Tail 
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